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Abstract 

The current fluctuations in normal metal / d-wave superconductor junctions 
are studied for various orientation of the crystal by taking account of the spa- 
tial variation of the pair potentials. Not only the zero-energy Andreev bound 
states (ZES) but also the non-zero energy Andreev bound states influence 
on the properties of differential shot noise. At the tunneling limit, the noise 
power to current ratio at zero voltage becomes 0, once the ZES are formed at 
the interface. Under the presence of a subdominant s-wave component at the 
interface which breaks time-reversal symmetry, the ratio becomes 4e. 
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The origin of shot noise is the current fluctuations in transport due to the discreteness 
of the charge carrier. Shot noise measurements provide important information on conduc- 
tion processes which can not be obtained from the usual conductance measurements^. In 
the last few years, several novel features peculiar to shot noise in mesoscopic systems have 
been revealed. In particular, the shot noise in normal metal-superconducting junction!! - ! 
and superconductor / insulator / superconductor junctional have been intensively stud- 
ied. It has been shown, through these works, that the Andreev reflection and the charge 
transport by the Cooper pairs have significant influence on the transport fluctuation at low 
voltages. However, most of previous theories are constructed on the conventional s-wave 
superconductors and the theory for <i-wave superconductors has not been presented. 

On the other hand, extensive experimental and theoretical investigations have revealed 
that the pair potentials of high-T c superconductors are d x 2_ y 2-wave symmetry.§ One of the 
essential differences of d^-j^-wave superconductors from conventional s-wave superconduc- 
tors is that the phase of the pair potential strongly depends on the wave vector. For example, 
the appearance of zero-bias conductance peak (ZBCP) in tunneling conductance at a (110) 
surface of d x 2_ y 2-wave superconductors reflects the sign change of effective pair potential 
through the reflection of quasiparticle at the surface.00 Orientational and material de- 
pendences of ZBCP of high-T c superconductors have been experimentally studied in several 
groups, and the consistency between theory and experiments has been checked in details&E. 

At this stage, it is an interesting problem to clarify what is expected in the shot noise in 
normal metal / insulator /d x 2 _ y 2 -wave superconductor (n/I/d) junction under the presence 
of the ZBCP. Recently, Zhu and Ting presented a theory of the shot noise in n/I/d junction^. 
They found a remarkable feature for n/I/d junction : when the angle between the normal 
to the interface a is ±7r/4, the noise-to-current ratio is zero at zero-bias voltage and quickly 
reaches a classical Shottky value 2e at finite voltage. This feature is completely discrepant 
from that for conventional s-wave superconductor (n/I/s) junctions where the ratio is 4e 
at zero voltage and 2e at finite voltage. This anomalous behavior is responsible for the 
formation of the zero energy Andreev bound states (ZES) at the interface of the <i-wave 



superconductor. Although Zhu and Ting theory (ZT theory) clarified important aspects of 
the shot noise under the presence of ZES, there still remain several unresolved problems. 

One is the detailed orientational dependence of the shot noise. In ZT theory, the vanish- 
ment of the noise-to-current ratio R(eV) at zero voltage (eV = 0) is shown only for a = tt/A 
(0 < a < vr/4). Since this property is related to the existence of ZBCP in tunneling con- 
ductance, we must check the value of -R(O) for < a < tt/4 where ZBCP appears in the 
tunneling conductance. In this paper, we will show that in the tunneling limit, i.e., low 
transparency limit, -R(O) vanishes for a^fl and becomes 4e only for a = where no ZES 
are expected. Moreover, it is shown that -R(O) is classified into three values, 0, 2e and 4e, 
corresponding to the region of the Fermi surface contributing to the ZES. 

The other is the influence of the spatial dependence of the pair potential on the shot 
noise. It is known that when the ZES are formed at the interface of <i-wave superconductor 
the pair potential is suppressed near the interfacelllfil. Consequently, not only the ZES but 
also the non-zero energy Andreev bound states (NZES) are formedfi The influences of the 
NZES on the shot noise are clarified. We further study the situation where a subdominant 
s-wave component which breaks the time reversal symmetry is induced near the interface 
of (i-wave superconductoiiiHil. The subdominant s-wave component influences significantly 
on R(eV). 

The model examined here is a two-dimensional n/I/d junction within the quasiclassical 
formalisml!@ where the pair potential has a spatial depends 



lcncc 



0, (x < 0) 

A(x,9) = { (1) 
A R (x,6), (x>0) 

Here 9 is the angle of quasiparticle trajectory measured from the x axis. If we apply this 
formula to d-wave superconductors including a subdominant s-wave component near the 
interface, A R (x,9) is decomposed into 

A R (x, 9) = A d {x) cos[2(# - a)} + A s (x) (2) 

where a denotes the angle between the normal to the interface and the x axis of the crystal. 



The insulator located between the normal metal and the superconductor is modeled by a 5 
function. The magnitude of the 5-function denoted as H determines the transparency of the 
junction a^, with = cos 2 9/[Z 2 + cos 2 #] and Z = mH /h 2 kp. The effective mass m and 
Fermi momentum /jp are assumed to be constant throughout the junction. The noise power 
to current ratio R(eV), the differential shot noise Sj>(eV), and the tunneling conductance 
as(eV), are given by@0 

Jo eV dES T (E) 
J dEa s {E) 

S T (eV) = - [ W 2 d6S(eV, 6) cos#, a s (eV) = - d6a s {eV, 6) cos# (4) 

2 J-7r/2 2 J-tt/2 

5(eV, 0) = ^-[Ra(l - Ra) + Rb(l - Rb) + 2RaRb] (5) 

Gs{E,e) = ^-{l + R a -R b ) (6) 



The magnitude of the Andreev and normal reflection R a and Rb are given b 

R = ^n\Vr,+ (°^)\ 2 (7) 

a |i + k-i)^, + (o,%b,-(o,#)I 2 ' 



and 
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R = (1 - Qjy) |1 -rj R ,+ ( Q ^)VR,-( Q ^j\ /g\ 

6 |i + K - 1)^(0, ^-(o^)! 2 " 

As a reference, we also calculate the tunneling conductance normalized by that in the normal 

state, 

aT(eV) = jTjW , (9) 

v 7 1 r n / 2 jo a v 7 

2 J-n/2 dv&N COSV 

Note that the differential shot noise and conductance spectrum are expressed only by 
t}r,±(x,9) just at the boundary {x = 0) where i] R ±(x,9) obeys the following equations, 
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—Vr,+ {x,u) = — 

dx invpcosu 



-&r(x, e+)v%,+(x, 6) - A R (x, 6+) + 2Er iRt+ (x, 9)\ , (10) 
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-A* R {x,e.)r J l_{x,e)-A R (x,e-)+2Eri R> -{x,e) , (11) 



dx ' ' ihvF cos 9 
with vp = kp/m, 9 + = 9 and 9_ = n — 9, and A R (x,9 + ) [A R (x, 6L)] is the effective 
pair potential felt by an electron [a hole] like quasiparticle. The quasiparticle energy E is 
measured from the Fermi energy. 

The spatial dependence of the pair potentials are determined by the following equations^ 



A s (x) = 9s k B Tj2^- d9'{[g R (9',x)] 12 - [g+{9' , x)] l2 } (12) 



A d (x) = g d k B Tj2 ^ I*'* d6'cos[2(6' - a)]{[g R (9',x)] l2 - \g^9',x)] 12 } (13) 



lim AJx) = 0, lim AJx) = A (14) 
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with dimensionless inter-electron potential of the s-wave g s and (i-wave g^, respectively. The 
quasiclassical Green's function g R (9,x) obeysS 

g R (9, x) = U R (9, x, 0)g R (9, 0)U R \9, x, 0) (15) 



d 

ihv Fx —U R (9,x, 0) 



t iuj n A R (x, 9 + ) ^ 



U R (9,x,0), (16) 



y -A* R (x, 9 + ) -iu n j 

with u n = l^ksTin + 1/2) and U R (9, 0,0) = 1. In the actual numerical calculations, 
f] R ±(x,9) is calculated from Eqs. ([K]) to (|TTD. Since g R (9,0) is expressed by rj R ±(9, 0), 
g R {9,x) is obtained using Eqs. (|15"D to (|T6|). Subsequently, the spatial dependence of the 
pair potentials A^(x) and A s (x) are calculated by Eqs. (|12|) to fli~4]). To get self-consistently 
determined pair potential, this process is repeated until enough convergence is obtained. 

First let us consider the case where A s (x) is not present. It is known for q ^0, apieV) 
has a zero bias conductance peak (ZBCP) with small magnitude of a It is expected 

that this property reflects on Sx(eV) for a ^ 0. In Fig. 1, SpieV) is plotted for Z = 5 with 
a = 7r/6 [see curve a in Fig. 1(a)]. As a reference, similar calculations are performed based 
on the non self-consistent pair potential where the spatial dependence of the pair potential 



is chosen as Aa(x) = A [see curve b in Fig. 1(a)]). As seen from curves a and b, St{0) has a 
peak around zero voltage originating from the ZBCP in or(0) [see Fig. 1(b)]. As compared 
to curves a to b, both the height and width of the peak around zero voltage of curve a are 
small as compared to those in b, since the degree of resonance at zero voltage is weakened 
due to the reduction of A d (x) at the interface, Besides this property, curves a in and Sy(eV) 
(<jt{g-V)) have second peak around eV ~ 0.45Ao due to the formation of NZES which can 
not be expected in curves b. 

In Fig. 2, we plot R(eV) for n/I/d junction for sufficient low transparency case, e.g., 
Z = 5. In the case of a = 0, no ZES are formed at the interface, then the resulting R(eV) is 
4e at zero voltage and is 2e at higher voltage (curve a in Fig. 2). When the value of a deviates 
from 0, R(0) is zero and it quickly reaches a classical Shottky value 2e, which is consistent 
with other reportS (curve b and c). This feature is explained based on the 9 dependence of 
a s (0,9) and S(0,6) as follows: in general, a s (0,9) = ^R a and S(0,6) = ^-RJX ~ R*) 
are satisfied. At the tunneling limit, R a is suppressed with the increase of Z unless the ZES 
are formed at the interface. When ZES are formed, R a = 1 is satisfied independent of Z. 
Such a situation is realized for it/A — a <\ 9 |< tt/A + a. For a = 0, no ZES are formed 
and both the magnitude of Sr(0) and crs(0) are reduced with the increase of Z, while R(0) 
remains to be constant. On the other hand, for a ^ the magnitude of St{0) is reduced 
while e>"s(0) remains finite with the increase of Z, thus -R(O) becomes zero. The important 
point is that the origin of the vanishment of -R(O) is responsible for the presence of ZES. It 
is a universal property excepted for unconventional superconductor junctions where finite 
region of the Fermi surface contributing on the formation of the ZES. 

The critical situation is realized in n/I/p x + ip y junction where ZES are formed only by 
the quasiparticles injected perpendicular to the interface^. The spatial dependence of the 
pair potential A R (x, 9) is given by 

A R (x,9) = A pl (x) cos9 + iA p2 (x) sin9 (17) 

with 
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lim A„i(x) = A , lim A p2 (x) = A . (18) 

R(eV) is calculated based on the self-consistently determined pair potentials, A pl (x) and 
A P 2(x). As shown in curve d in Fig. 2, R(eV) is nearly 2e independent of bias voltages. 
The nature that -R(O) is neither 4e nor has never been predicted in previous theorieslE. 

Under the presence of the ZBCP, since the quasiparticle density of states at the zero en- 
ergy near the interface are enhanced, a subdominant s-wave component of the pair potential 
A s (x) can be induced near the interface, when a finite s-wave pairing interaction strength 
exists, even though the bulk symmetry remains pure <i-wave§iHil. Since the phase difference 
of the d-wave and s-wave components is not a multiple of 7T, the mixed state breaks the 
time-reversal symmetry§!I!. In such a case, the ZBCP splits into two and the amplitude of 
the splitting depends on the magnitude of the induced s-wave component. The correspond- 
ing R(eV) is plotted in Fig. 3(a) with a = 7r/4 and Z = 5 where the transition temperature 
of s-wave component is chosen as T s = 0.15T^ (dotted line) and T s = 0.3T^ (solid line). 
The induced s-wave component near the interface influences crucially on the R(eV) at low 
voltages. The remarkable feature is that -R(O) recovers to be 4e as in the cases of n/I/s 
junctions. It is because that with the inducement of s-wave component which breaks the 
time reversal symmetry, the position of the Andreev bound state shifts to eV ~ A s (0), 
where ar(eV) has a peak [see Fig. 3(b)]. Consequently, ZBCP shifts to this voltage and 
R(eV) has a dip structure. Although -R(O) = 4e is satisfied, the overall feature of R(eV) is 
completely different from that in n/I/s junction, since s-wave component is induced only 
near the interface in the present case. 

In this paper, the current fluctuations in normal metal / c?-wave superconductor junctions 
are studied for various orientation of the junction by taking account of the spatial variation 
of the pair potentials. Not only the zero energy Andreev bound states but also the non-zero 
energy Andreev bound states show up in the line shape of St(gV). At the tunneling limit, 
we found universal property of -R(O) for two dimensional superconductors. -R(O) is 0, 2e and 
4e, corresponding to three cases where the region of the Fermi surface contributing to the 
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ZES is i)finite region, ii)point and iii)none, respectively. The present property gives useful 
information on the identification for the symmetry of the unconventional superconductors. 
We hope the measurement of shot noise in unconventional superconductors will be performed 
near future. 

It is known that ZES also influences significantly on Josephson current in d-wave super- 
conductor / insulator / <i-wave superconductor (d/I/d) junctionsil'il. It is an interesting 
and future problem to study the shot noise in d/I/d junctions. 
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FIG.l 

FIG. 1. a T (eV) [Fig.l(a)] and S T {eV) [Fig.l(b)] for n/I/d junction with Z = 5 and a = vr/6. 
a: self-consistent calculation, b: non self-consistent calculation. 
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FIG.2 

FIG. 2. R(eV) is plotted with Z = 5 a: a = 0b: a = n/12 c: a = ir/4 for n/I/d junctions, d: 
Similar plot for n/I/p junction. 
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FIG. 3 

FIG. 3. R(eV) is plotted for Z = 5 and a = ir/A where s-wave component is induced at the 
interface [Fig. 3(a)]. a: T s = 0.15^ and b: T s = 0.3^. The corresponding ax(eV) is plotted in 
[Fig. 3(b)]. 
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